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Abstract 

We calculate the cross section of beauty production in ep deep inelastic scattering at 
HERA collider in the framework of the fcT^-factorization approach. The unintegrated gluon 
distributions in a proton are obtained from the full CCFM, from unified BFKL-DGLAP 
evolution equations as well as from the Kimber-Martin-Ryskin prescription. We investigate 
different production rates and study the 6-quark contribution to the inclusive proton struc- 
ture function F2(x, Q^) at small x and at moderate and high values of Q^. Our theoretical 
results are compared with the recent experimental data taken by the HI and ZEUS collab- 
orations. We demonstrate the importance of leading Inl/x contributions in description of 
the HERA data. 



1 Introduction 

The beauty production at high energies is a subject of intensive study from both theo- 
retical and experimental points of view [1-9]. First measurements [1] of the 6-quark cross 
sections at HERA were significantly higher than the QCD predictions calculated at next- 
to-leading order (NLO) approximation. Similar observations were made in hadron-hadron 
collisions at Tevatron [2] and also in photon-photon interactions at LEP2 [3]. In last case, 
the theoretical NLO QCD predictions are more than three standard deviations below the 
experimental data. At Tevatron, recent analisys indicates that the overall description of the 
data can be improved [10] by adopting the non-perturbative fragmentation function of the 
6-quark into the i?-meson: an appropriate treatment of the 6-quark fragmentation properties 
considerably reduces the disagreement between measured beauty cross section and the results 



1 



of corresponding NLO QCD calculations. Also latest measurements [4, 5, 9] of the beauty 
photoproduction at HERA are in a reasonable agreement with the NLO QCD predictions 
or somewhat higher. Some disagreement is observed mainly at small decay muon and/or 
associated jet transverse momenta [4, 5, 9]. But the large excess of the first measurements 
over NLO QCD, reported by the HI collaboration [1], is not confirmed. 

Recently there have been become available experimental data [6-9] on the 6-quark produc- 
tion in deep inelastic scattering at HERA which taken by the HI and ZEUS collaborations. 
The first measurements [6, 7] of the beauty contribution to the inclusive proton structure 
function F2{x, Q^) have been presented for small values of the Bjorken scaling variable x, 
namely 2 • 10~^ < x < 5 ■ 10~^, and for moderate and high values of the photon virtuality 
Q^, namely 12 < < 700 GcV^. Also process e + p^e' + 6 + 6 + X->e' + jet + /i + X' 
has been measured [8, 9] in the small x region with at least one jet and a decay muon in the 
final state and still was not described in the framework of QCD theory. Such processes are 
dominated by the photon-gluon fusion subprocess 7* + — > 6 + 6 and therefore sensitive to 
the gluon density in a proton xg{x, fi'^). It was claimed [8, 9] that the NLO QCD calculations 
have some difficulties in description of the recent HERA data. The predictions at low values 
of Q^, Bjorken x, muon transverse momentum and high values of jet transverse energy and 
muon pseudo-rapidity is about two standard deviation below the data. 

In the present paper to analyze the recent HI and ZEUS data [6-9] we use the so-called 
/cT-factorization [11, 12] (or semi-hard [13, 14]) approach of QCD which has been applied 
earlier, in particular, in description of the charm and beauty production at HERA [15-21], 
Tevatron [22-28] and LEP2 [19, 29, 30] colliders. The /cr-factorization approach is based 
on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [31] or Ciafaloni-Catani-Fiorani-Marchesini 
(CCFM) [32] gluon evolution which are valid at small x since here large logarithmic terms 
proportional to Inl/x are summed up to all orders of perturbation theory (in the leading 
logarithmic approximation). It is in contrast with the popular Dokshitzer-Gribov-Lipatov- 
Altarclli-Parizi (DGLAP) [33] strategy where only large logarithmic terms proportional to 
In/x^ are taken into account. The basic dynamical quantity of the /cy-factorization approach 
is the so-called unintegrated (k^-dependent) gluon distribution A{x, kf., /i^) which deter- 
mines the probability to find a gluon carrying the longitudinal momentum fraction x and 
the transverse momentum k^ at the probing scale /i^. The unintegrated gluon distribution 
can be obtained from the analytical or numerical solution of the BFKL or CCFM evolution 
equations. Similar to DGLAP, to calculate the cross sections of any physical process the un- 
integrated gluon density A{x, k^, //^) has to be convoluted [11-14] with the relevant partonic 
cross section a. But as the virtualities of the propagating gluons are no longer ordered, the 
partonic cross section has to be taken off mass shell (kr-dependent). It is in clear contrast 
with the DGLAP scheme (so-called collinear factorization). Since gluons in initial state are 
not on-shell and are characterized by virtual masses (proportional to their transverse mo- 
mentum), it also assumes a modification of their polarization density matrix [13, 14]. In 
particular, the polarization vector of a gluon is no longer purely transversal, but acquires 
an admixture of longitudinal and time-like components. Other important properties of the 
/^T-factorization formalism are the additional contribution to the cross sections due to the 
integration over the k|, region above jj,^ and the broadening of the transverse momentum 
distributions due to extra transverse momentum of the colliding partons. 

As it was noted already, some applications of the /cT-factorization approach supplemented 
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with the BFKL and CCFM evolution to the heavy (charm and beauty) quark production 
at high energies are widely discussed in the literature [15-30] (see also review [34, 35]). It 
was shown [24-28] that the beauty cross section at Tevatron can be consistently described 
in the framework of this approach. However, a substantial discrepancy between theory and 
experiment is still found [19, 29, 30] for the &-quark production in 77 collisions at LEP2, 
not being cured by the fcr-factorization^. At HERA, the inclusive beauty photoproduction 
has been investigated [16, 17, 19, 21, 25]. In [17, 19, 25] comparisons with the first HI 
measurements [1] have been done. In [17, 25] the Monte-Carlo generator Cascade [36] has 
been used to predict the cross section of the 6-quark and dijet associated photoproduction. 
However, calculations [17, 19, 25] deal with the total cross sections only. In our previous 
paper [21] the total and differential cross sections of beauty photoproduction (both inclusive 
and associated with hadronic jets) have been considered and comparisons with the recent 
HI and ZEUS measurements [1, 4, 5, 9] have been made. It was demonstrated [21] that 
the /cT-factorization approach supplemented with the CCFM or BFKL-DGLAP evolved un- 
integrated gluon distributions [28, 37] reproduces well the numerous HERA data [1, 4, 5, 
9]. 

In the present paper we will study the beauty production in ep deep inelastic scatter- 
ing at HERA. We investigate a number of different production rates (in particular, the 
transverse momentum and pseudo-rapidity distributions of muons which originate from the 
semi-leptonic decays of 6-quarks). Our study is based on leading-order (LO) off-mass shell 
matrix elements for the photon-gluon fusion subprocess e + g* ^ e' + b + b. Particularly 
we discuss the photoproduction limit (Q^ 0) of our derivation. Also we investigate the 
beauty contribution to the inclusive proton structure function F2{x,Q^). In the numerical 
analysis we test the unintegrated gluon distributions which were obtained [28, 37] from the 
full CCFM, unified BFKL-DGLAP evolution equations and from the conventional parton 
densities (using the Kimbcr-Martin-Ryskin prescription [38]). We attempted a systematic 
comparison of model predictions with the recent experimental data [6-9] taken by the HI 
and ZEUS collaborations. One of purposes of this paper is to investigate the specific kx- 
factorization effects in the 6-quark leptoproduction at HERA. 

The outline of our paper is following. In Section 2 we recall the basic formulas of the 
^T-factorization approach with a brief review of calculation steps. In Section 3 we present 
the numerical results of our calculations and a discussion. Finally, in Section 4, we give some 
conclusions. The compact analytic expressions for the off-mass shell matrix elements of the 
photon-gluon fusion subprocess e + g* ^ e' + b + b are given in Appendix. These formulas 
may be useful for the subsequent applications. 

2 Calculation details 

In this section we present our analytic results for the cross section of e+p — * e' -\-b + b + X 
in DIS. We work at leading-order /cr-factorization approach of QCD. We start by defining 
the kinematics. 

^Some discussions of this problem may be found in [19, 30]. 
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2.1 Kinematics 



We denote the four-momenta of the incoming electron and proton and the outgoing 
electron, beauty quark and anti-quark by pe, Pp, p'^, Pb and pi, respectively. The off-shell 
gluon and virtual photon have four-momenta k and q, and it is customary to define 

= 4 = -4 < 0, = (pe - p'f ^q^T- -Q^ < 0, (1) 

where kx and qx are the transverse four-momenta of the corresponding particles. Choosing 
a suitable coordinate system in the ep center-of-mass frame, we have 

Pe = %/5/2 (1,0, 0,-1), Pp = x/5/2(l,0,0,l), (2) 

where is the total energy of the process under consideration and we neglect the masses of 
the incoming electron and proton. The standard deep inelastic variables x and y are defined 
as usual: 

Q' {Pp ■ q) Q' .ON 

2(Pp-g) KPe-Pp) xs 

The variable y measures the relative electron energy loss in the proton rest frame. Prom the 
conservation law we can easily obtain the following condition: 

kr + Qt = P6T + VbT- (4) 

2.2 Cross section for deep inelastic beauty production 

According to the /cr-factorization theorem, the cross section of deep inelastic beauty 
production e -\- p ^ e' -\- h -\-h -\- X can be written as a convolution 

dd) — 
— A{x, 14, ^l'^)dkl,—dc^{e + g* ^ e' + h + h), (5) 
X 2tt 

where A{x, k"^, 11^) is the unintegrated gluon distribution in a proton, a{e + g* — >• e' -1-6-1-5) is 
the cross section of partonic subprocess and (f) is the azimuthal angle of initial virtual gluon. 
Decomposing the cross section &{e -\- g* ^ e' -\- h -\-h) into a leptonic and a hadronic part, 
we can write it as 

da{e + g*^e' + b + b)^ _l_^L^-if^^ rf$(3)(p^ + k,p'^,Pb,Pi), (6) 

where e is the electron charge magnitude and L^^'^ and H^^'^ are the leptonic and hadronic 
tensors. In general case the Lorentz- invariant element . . . ,p„) of n-body phase 

space is given by 

...,p^) = (2vr^) 5W - ±p^ n (^|%o- (7) 
Integrating over the azimuthal angle of the outgoing electron, we can simplify (6) to become 
da(e + g*^e' + b + b) = ^^^^^"^ H,,J^d^^'\^ + k,Pb,Pi), (8) 



where ol = e^/(47r) is Sommerfeld's fine structure constant. For the leptonic tensor L'^'^ we 
use the following expression 



^ ^ + (^-y)\.^ _ (9) 

y y 

where 

The ej.'' and e'^ refer to transverse and longitudinal virtual photon polarization, as indicated 
by their subscripts. It is easily to see that q^e^y — q^e'^ — 0, e'^j^ — —2 and eJ^J;^ = — 1. 
Furthermore, 

e^^-^^^ + ef = + (11) 

i.e. e'^'^ is the polarization tensor of an unpolarized spin-one boson having mass q^. From (5) 
— (10) one can obtain the following formula for the cross section of deep inelastic beauty 
production in the /c^-factorization approach: 

,(,+^^e' + 6 + 5 + X) = / k^, ^^') X ^^^^ 



X 



2 Ji 2 r , d4)d(t>bd(t)i 



dpf^j^dkj^dQ dybdyi 



27r 27r 27r ' 



where y^, m and (ph and (pi are the rapidities and azimuthal angles of the produced beauty 
quark and anti-quark, respectively. The evaluation of functions TCk"^, Q'^) and L{\i\, Q'^) 
has been done analytically using the Mathematica 5 program. The compact expressions 
for these functions are listed in Appendix. It is important that the functions T(k|,, Q^) 
and L(k|., Q"^) depend on the virtual gluon non-zero transverse momentum k^. Note that 
if we average (12) over k^ and take the limit k|, 0, then we obtain well-known formula 
corresponding to the usual LO QCD calculations. 

It is interesting to study the photoproduction limit of (12) by taking the limit — > 0. 
This provides us with a powerful check for our formulas by relating them to well-known 
results. So, the cross section of the partonic process 7 -|- gi* — > 6 -|- 6 reads 



d^^-'\q + k,p,,pi), (13) 



where s — {q-\- k) . Comparing (8) and (13), one can obtain the well-known relation 



^da{e + g* ^ e' + h + b) a I + {1 - yf 



jgo^ \lyaQ^ ^ = ^^^-(7 + /-^ + ^. (14) 

The contribution of 6-quarks to the deep inelastic proton structure function ^2(2;, Q"^) can be 
calculated according to convolution (5) also. The relevant coefficient function is described 
by the quark box diagram and has been presented in our previous paper [40] . 
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The multidimensional integration in (12) has been performed by means of the Monte 
Carlo technique, using the routine Vegas [41]. The full C++ code is available from the 
authors on request^. 

3 Numerical results 

We now are in a position to present our numerical results. First we describe our theoret- 
ical input and the kinematical conditions. 

3.1 Theoretical uncertainties 

There are several parameters which determined the normalization factor of the cross 
section (12): the beauty mass m^, the factorization and normalisation scales iip and and 
the unintegratcd gluon distributions in a proton A{x^ k|,,yU^). 

Concerning the unintegrated gluon densities in a proton, in the numerical calculations 
we used five different sets of them, namely the J2003 (set 1 — 3) [28], KMS [37] and 
KMR [38]. All these distributions are widely discussed in the literature (see, for example, 
review [34, 35] for more information). Here we only shortly discuss their characteristic 
properties. First, three sets of the J2003 gluon density have been obtained [28] from the 
numerical solution of the full CCFM equation. The input parameters were fitted to describe 
the proton structure function F2{x, Q"^). Note that the J2003 set 1 and J2003 set 3 densities 
contain only singular terms in the CCFM splitting function Pgg{z). The J2003 set 2 gluon 
density takes into account the additional non-singlular terms'^. These distributions have 
been applied in the analysis of the forward jet production at HERA and charm and bottom 
production at Tevatron [28] (in the framework of Monte-Carlo generator Cascade [36]) and 
have been used also in our calculations [20, 21]. 

Another set (the KMS) [37] was obtained from a unified BFKL-DGLAP description of 
F2{x,Q'^) data and includes the so-called consistency constraint [42]. The consistency con- 
straint introduces a large correction to the LO BFKL equation. It was argued [43] that about 
70% of the full NLO corrections to the BFKL exponent A are effectively included in this con- 
straint. The KMS gluon density is successful in description of the beauty hadroproduction 
at Tevatron [24, 26] and photopro duct ion at HERA [21]. 

The last, fifth unintegrated gluon distribution A{x, k|., /i^) used here (the so-called KMR 
distribution) is the one which was originally proposed in [38]. The KMR approach is the 
formalism to construct unintegrated gluon distribution from the known conventional parton 
(quark and gluon) densities. It accounts for the angular-ordering (which comes from the 
coherence effects in gluon emission) as well as the main part of the collinear higher-order 
QCD corrections. The key observation here is that the fi dependence of the unintegrated 
parton distribution enters at the last step of the evolution, and therefore single scale evolution 
equations (DGLAP or unified BFKL-DGLAP) can be used up to this step. Also it was 
shown [38] that the unintegrated distributions obtained via unified BFKL-DGLAP evolution 
are rather similar to those based on the pure DGLAP equations. It is because the imposition 

^ lipatov® theory, sinp.msu.ru 
^See Ref. [28] for more details. 
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of the angular ordering constraint is more important [38] than inchiding the BFKL effects. 
Based on this point, in our further calculations we use much more simpler DGLAP equation 
up to the last evolution step^. Note that the KMR parton densities in a proton were used, 
in particular, to describe the prompt photon photoproduction at HERA [45] and prompt 
photon hadroproduction Tevatron [46, 47]. 

Also the significant theoretical uncertainties in our results connect with the choice of 
the factorization and renormalization scales. First of them is related to the evolution of 
the gluon distributions, the other is responsible for the strong coupling constant ^^(/xfj). 
The optimal values of these scales are such that the contribution of higher orders in the 
perturbative expansion is minimal. As it often done for beauty production, wc choose the 
renormalization and factorization scales to be equal: = = ^ = ^ ml + (py ) , where 
(p^) is set to the average p^, of the beauty quark and antiquark. But in the case of the 
KMS gluon distribution we used special choice it was originally proposed in [37] . 

Note that in the present paper we concentrate mostly on the non-coUinear gluon evolution 
in the proton and do not study the scale dependence of our results. To completeness, we 
take the 6-quark mass = 4.75 GcV and use LO formula for the coupling constant ^^(/x^) 
with Uf — 4 active quark flavours at Aqcd = 200 MeV, such that as{M'^) — 0.1232. 

3.2 Associated beauty and jet production 

The recent experimental data [8, 9] for the associated beauty and hadronic jet Icpto- 
production at HERA comes from both the HI and ZEUS collaborations. The total and 
differential cross sections as a function of the photon virtuality Q^, Bjorken scaling variable 
X, muon transverse momentum and pseudo-rapidity rj'^ and jet transverse momentum 
Pj^ have been determined. The ZEUS data [8] refer to the kinematical region^ defined by 

> 2 GeV^ with at least one hadron-level jet (in the Breit frame) with p^^*^"^*^'* > 6 GeV 
and —2 < rf'^^ < 2.5 and with muon which fulfill the following conditions: —0.9 < 1]^^ < 1.3 
and pj. > 2 GeV or —1.6 < rj^ < —0.9 and > 2 GeV. The fraction y of the electron energy 
transferred to the photon is restricted to the range 0.05 < y < 0.7. Note that the Breit frame 
is defined by the usual condition q+2xpp — 0. In this frame, a space- like photon and proton 
collide head-on and any final-state particle with a high transverse momentum is produced by 
a hard QCD interaction. The more recent HI data [9] refer to the kinematical region defined 
by 2 < g2 < 100 GeV2, Q.K y < 0.7, > 2.5 GeV, -0.75 < 7]^" < 1.15, y^*^'"* > 6 GeV 
and < 2.5. To produce muons from 6-quarks in our theoretical calculations, we first 
convert 6-quarks into S-hadrons using the Peterson fragmentation function [48] and then 
simulate their semileptonic decay according to the standard elcctrowcak theory. Of course, 
the muon transverse momenta spectra are sensitive to the fragmentation functions. However, 
this dependence is expected to be small as compared with the uncertainties coming from the 
unintegrated gluon densities in a proton. Our default set of the fragmentation parameter is 
66 = 0.0035. 

The basic photon-gluon fusion subprocess '-/* + g* ^ b + b give rise to two high-energy 
6-quarks, which can further evolve into hadron jets. In our calculations we assumed that the 

^We have used the standard GRV (LO) parametrizations [44] of the coUincar quark and ghion densities. 
^Here and in the following all kinematic quantities are given in the laboratory frame where positive OZ 
axis direction is given by the proton beam. 
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produced quarks (with their known kinematical parameters) are taken to play the role of the 
final jets. These two quarks are accompanied by a number of gluons radiated in the course 
of the gluon evolution. As it has been noted in [15], on the average the gluon transverse 
momentum decreases from the hard interaction block towards the proton. We assume that 
the gluon emitted in the last evolution step and having the four-momenta k' compensates 
the whole transverse momentum of the gluon participating in the hard subprocess, i.e. kV — 
— ky. All the other emitted gluons are collected together in the proton remnant, which is 
assumed^ to carry only a negligible transverse momentum compared to k'y This gluon 
gives rise to a final hadron jet with p^* = |kV| in addition to the jet produced in the hard 
subprocess. From these three hadron jets we choose the one carrying the largest transverse 
momentum (in the Breit frame), and then compute the beauty and associated jet production 
rates. 

The results of our calculations are shown in Figs. 1 — 10 in comparison to the HI and 
ZEUS experimental data [8, 9] fo the 6-quark and associated jet production. Sohd, dashed, 

dash-dotted, dotted and short dash-dotted curves correspond to the predictions obtained 
with the J2003 set 1 — 3, KMR and KMS unintegrated gluon densities, respectively. One 
can see that the overall agreement between our results (calculated using the J2003 and KMS 
gluon densities) and experimental data [8, 9] is a rather good. However, the measured cross 
section as a function of the muon transverse momentum pj. shows a slightly steeper behaviour 
than the theoretical predictions: the results of our calculations tends to underestimate the 
data at low pj. (see Figs. 1 and 2). But in general these predictions still agree with the 
HI and ZEUS data within the experimental uncertainties. Note also that the measured 
differential cross sections da/drj'^ in Figs. 3 and 4 exhibit a rise towards the forward muon 
pseudo-rapidity region, which is not reproduced [8, 9] by the coUinear NLO calculations. 
At the same time the shape and the normalization of r/'^ distributions are well described by 
our calculations. The collinear NLO QCD underestimate also the data at low and low 
X values: it was claimed that in these kinematical regions the data are about two standard 
deviation higher [8, 9]. 

As it was already mentioned above, the absolute normalization of the predicted cross 
sections in the framework of fcr-factorization approach is depends on the unintegrated gluon 
distribution used. From Figs. 1 — 10 one can see that all three sets of the J2003 gluon 
density as well as the KMS one give rise to results which are rather close to each other. So, 
the difference in normalization between the KMS and J2003 predictions is rather small, is 
about 15% only. The similar effect we have found [21] in the case of beauty photoproduction. 
However, it is in the contrast with the D* meson and dijet associated photoproduction, which 
has been investigated in our previous paper [20]. It was demonstrated [20] a relative large 
enhancement of the cross sections calculated using the KMS gluon density. The possible 
explanation of this fact is that the large 6-quark mass (which provide a hard scale) makes 
predictions of the perturbation theory of QCD more applicable. Note also that the KMS 
gluon density provides a more hard transverse momentum distribution of the final muon (or 
jet) as compared with other unintegrated densities under consideration. Similar effect we 
have observed [21] in the case of beauty photoproduction. 

Concerning the KMR predictions, one can see that this unintegrated gluon distribution 

^Note that such assumption is also used in the KMR formahsm. 
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Source 


a{e + p ^ e' + jet + fi + X) [pb] 


ZEUS measurement [8] 
NLO QCD (hvqdis [50]) 
Rapgap [48] 
Cascade [36] 

jzuuo sei 1 
J2003 set 2 
J2003 set 3 
KMR 
KMS 


40.9 ±5.7 (stat.)+4 4 (syst.) 

20.6t|^ ■ 

14.0 

28.0 
'^p; 97 

33.47 
36.75 
22.11 
38.52 



Table 1: The total cross section of beauty and associated jet leptoproduction obtained in 
the kinematic range > 2 GeV^, 0.05 <y < 0.7, > 6 GeV, -2 < r]^^^ < 2.5 and 

> 2 GeV, -0.9 <r]^' < 1.3 or > 2 GeV, -1.6 < r]^' < -0.9. 

gives results which lie below the data and which are very similar to the coUinear NLO QCD. 
Such observation coincides with the ones [20, 21]. This fact confirms the assumption which 
was made in [45] that the KMR formalism results in some underestimation of the predicted 
cross sections. Such underestimation can be explained by the fact that leading logarithmic 
terms proportional to \nl/x are not included into the KMR approach. 

Now we turn to the total cross section of 6-quark and associated jet leptoproduction. In 
Table 1 and 2 we compare our theoretical results with the HI and ZEUS data [8, 9] obtained 
in relevant kinematical regions (defined above) . The predictions of Monte-Carlo generators 
Rapgap [49], Cascade [36] as well as NLO QCD calculations (hvqdis program) [50] are 
also shown for comparison. One can see that the collinear NLO QCD predictions is about 2.5 
standard deviation lower than the ZEUS data and is about 1.8 standard deviation lower then 
the HI data. At the same time, our predictions obtained using the J2003 and KMS gluon 
densities are significantly higher and agree well with the both HI and ZEUS data within 
the experimental uncertainties. The KMR unintegrated gluon distribution again gives the 
results which are below the data and which are very close to NLO QCD ones. Note that the 
Monte-Carlo generators Rapgap and Cascade also predict a lower cross section than that 
measured in the data. 

In general, we can conclude that the cross sections of deep inelastic beauty and associated 
jet production calculated in the A;7--factorization formalism (supplemented with the CCFM 
or unified BFKL-DGLAP evolution) are larger by 30 — 40% than ones calculated at NLO 
level of collinear QCD. This enhancement comes, in particular, from the non-zero transverse 
momentum of the incoming off-shell gluons and from taken into account the leading In 1 /x 
terms. Our results for the total and differential cross sections are in a better agreement (both 
in normalization and shape) with the HI and ZEUS data than the NLO QCD predictions. 

3.3 Beauty contribution to the proton SF F2{x, Q'^) 

Now we will concentrate on the 6-quark contribution to the inclusive proton structure 
function F2{x,Q^). We will use the master formulas which were obtained in our previous 
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Source 


a{e + p ^ e' + jet + fi + X) [pb] 


HI measurement [9] 
NLO QCD (hvqdis [50]) 
Rapgap [49] 
Cascade [36] 

jzuuo bei 1 
J2003 set 2 
J2003 set 3 
KMR 
KMS 


16.3 ± 2.0 (stat.) ± 2.3 (syst.) 
9.0l?j 
6.3 
9.8 

18.98 
20.80 
12.45 
22.61 



Table 2: The total cross section of beauty and associated jet leptoproduction obtained in 
the kinematic range 2 < < 100 GeV^, 0.1 < y < 0.7, > 2.5 GeV, -0.75 < ry^^ < 1.15, 

paper [40]. As it was mentioned above, the first experimental data [6, 7] on the structure 
function F2{x, Q'^) comes from the HI collaboration. These data refer to the kinematical 
region defined by 2 • 10"^ <x<b- IQ'^ and 12 < < 650 GeV^. 

Note that we change now the default set of parameters which we have used in the previous 
section. So, we set the renormalization and factorization scales /ir and /i^? to be equal to 
photon virtuality , as it was done earlier in analysis [5 1] of the charm contribution to the 
structure function ^2(2;, Q^) in the framework of /cT-factorization QCD approach. The similar 
choice have been used also in the analysis of longitudinal structure function Fl{x, Q^) [52]. 
Of course, in the case of the KMS gluon distribution we set = /J'^ — k|., as it was 
originally proposed in [37]. Other parameters have not been changed. 

In Fig. 11 we show the structure function ^2(0;, Q^) as a function of x for different values 
of in comparison to the recent HI data [6, 7]. One can see that the J2003 distributions 
reproduce well the experimental data for all values of Q^. The KMS gluon density demon- 
strates a perfect agreement with the data at moderate but slightly overestimate them at 

= 650 GeV^. It is interesting to note that the KMR density docs not contradict the ex- 
perimental data, too. However, this distribution predicts a more rapid rise of the calculated 
function F|(a;, Q^) with decreasing of x (in comparison to the J2003 and KMS densities). 
We can conclude that in the small x region {x < 10~^) the shape of function F2{x, Q^) 
predicted by the unintegrated gluon distributions under consideration is very different. In 
particular, the differences observed between the curves are due to the different behaviour of 
the corresponding unintegrated gluon distributions as a function of x and [45]. This fact 
shows the importance of a detail understanding of the non-coUinear parton evolution in a 
proton and the necessarity of better experimental constraints as well as further theoretical 
studies in this field. 

4 Conclusions 

We have calculated the deep inelastic beauty and associated jet production in electron- 



10 



proton collisions at HERA in the /cT^-factorization QCD approach. The total and several 
differential cross section (as a function of the photon virtuality Q^, Bjorken scaling variable 
X, decay muon transverse momentum pj, and pseudo-rapidity r]^ and hadronic jet transverse 
momentum pip') have been studied. Additionally we have investigated the 6-quark contribu- 
tion to the inclusive proton structure function F2{x, Q"^) at small x and at moderate and high 
Q^. In numerical analysis we have used the unintegrated gluon densities which are obtained 
from the full CCFM (J2003 set 1 — 3), from unified BFKL-DGLAP evolution equations 
(KMS) as well as from the Kimber-Martin-Ryskin prescription. Our investigations were 
based on the LO off-mass shell matrix elements for photon-gluon fusion subprocesses. 

We have shown that the A'^'-factorization approach supplemented with the CCFM or 
BFKL-DGLAP evolved unintegrated gluon distributions (the J2003 or KMS densities) re- 
produces well the numerous HERA data on beauty and associated jet production. At the 
same time we have obtained that the Kimber-Martin-Ryskin formahsm results in some un- 
derestimation of the cross sections. This shows the importance of a detail understanding of 
the non-coUinear parton evolution process. 
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6 Appendix 

Here we present the compact analytic expressions for the functions T(k|., Q^) and 
L(k|., Q^) which appear in (12). In the following, s, t and u are usual Mandelstam vari- 
ables for corresponding '-^* + g* ^ b + b subprocesses {s + t + u = 2w? — Q"^ — k|.) and m 
and Cfe is the mass and fractional electric charge of 6-quark. The exact expressions for the 
functions T(\<;^, Q"^) and L(k^, Q"^) can be presented as 

T(k|, = {An)Wa,{^,-)el ' f j,, ^ , (^-2) 

8(t — m^Y[u — m^y[s + -\- ky)^ 



FL{s,i,u,\4,Q'') 

8{i - m'^fiu - m2)2(s + + kf, )4 ' 



where 



Fris, t, u, k^, Q2) = -8(4k|g^(£ - + 2]4Q\-8m^ + + 4Q\t - uf- 

ii^u - 2pu'^ - 4iu^ + -a^ + 16m^{i + u)+ AQ'^{i - uf{i +u)- 8m'^{t^+ 
Aiu + u^) - Sm^{Q^{i - uf - 2iu{i + u)))+ 
(s + + k|)2(24m^^ + 8m^°(Q^ - 3(£ + -?/)) - 2m^(2Q^ + 3?+ 
22iu + Zv? + 10g2(t + u)) - iu{P + u^){2Q^ + 2Q^{i + u) + {i + 11)^)+ 



11 



4:QHu{3P + Aiu + Sm^)) + 4m^(2Q^(t + u) + Q^{3P + Uiu + 3u'^) + 
4(t^ + 6pM + 6tM^ + M^)) - m^(7f + 56Pxi + 90Pm^ + 56tM^ + 7u^+ 
6Q\i + uf + 2g2(f3 + i9f2^ + i9£^2 ^ ^3^^^ _ 2k^(8m^^ - 2Q^{i - uf- 

AQ\i - uf{i + u)+ iu{i + uYiP + u'') - AQ\i^ - 2Pu - 2Pu' - 2ie + u^) 
Q^(i^ - 5i% - At^ii^ - - 5iu'^ + u^) + 8m^°(2g2 - 3{i + u)) + 2m^(8Q^ 
12Q^{i + u) + 15{i + uf) - 4:m^{SQ^{t + u) + 5{i + uf+ 
g^(-5? + 2iu - 5u^)) + m^{7i^ + 32t^u + A2t^u'^ + 32iu^ + 7u^+ 
SQ'^iP + lOiu + u^) - 2(5^(5^ - I3pu - 13tu^ + 5u^)) + m^{-t^+ 
8Q%i - - 9i% - UPu^ - Up{? - 9tu;'^ - + 8Q\i^ - 5Pu- 
biu^ + u^) + AQ\i* - 5Pu - APv^ - 5iu^ + u'^))) + k|(32m^^- 
112m^2(£ + m) - 8m^°(4g'^ - 17P - 50iu ~ m?)- 
2iu{i + uf{t^ + u^) + Q2(P - u^f{P - Aiu + u^) + 
2Q\i'^ - 4Pu - 2Pu^ - Aiu^ + u^) + 2Q\P - 
APu^ - APu^ - 5iu^ + u^) - 4m^(4g6 + 19P- 
lAQ^{i - uf + 12lPu + 12liu'^ + 19u^ - 12Q^{i + u)) + Am^{5p+ 
72Pu + 126Pii^ + 72iu^ + ^u^ + SQ^{i +u)- 20Q^{i - uf{i + u)- 
2Q^{5P - 2iu + hu^)) - 2m^{P + A9Pu + ll%Pu^ + 118?# + 49M^+ 
+ 8g^(? + Aiu + u^) - 3Q'^{i - uf{7P + mil + 7u^)- 
2Q^{5P - 13Pu - 13iV + 5#)) + 2m'^{16QHu{i + u)+ 
2iu{i + uf{5p + Aiu + 5u^) - Q^{i - uf{5p + 3Pu+ 

3ie + 5^) - AQ^i" - 5pu - Ape - 5ie + «^)))), 

Fl{s, i, u, \4, Q^) = lQ{2\4Q^{i - uf + ]4{i - uf{2m^ + AQ^ + P + u^+ 
AQ^{i + u) - 2m^{AQ^ + t + m)) + 2{m^ - i){m^ - + + Wj>f{2m^^ 
m^{Q^ -i-u)+ iu{Q^ + i + u) - m^{2i{i + Q^{i + fi))) + \4{-8m^Q^+ 
2Q%i -uf + AQ\i - u)\i + ii) + a - u)\i + u) + 
Q2(3t^ - 6pM - 2Pu^ - 6iu^ + 3m^) + 8m^{-{i - uf + 2Q^{i + u))- 
2m^(-4(t - uf{i + u) + Q^{P + 22£m + u'^)) - 2m^{AQ^{i -uf+ 
2{i - fif{P + fj^) + Q^{3P - ni% - lltV + 3u^))) + 2\4{Am^^- 
4m^°(2Q^ + 3{i + u)) + m\-AQ^ + 17P + 26iu + 17u'^ + 12Q^{i + u))+ 

2m%AQ\i + u) - 5{i + uf - Q\P + 6iu + u^)) - iu{2Q\P + e)+ 
{i + uf{P - 3iu + u^) + Q'^{3P + Pu + iu^ + 3#)) + 2m^{Q^{i + uf+ 
iu{P - 7Pu - 7iu^ + + g^(£^ + hPu + biu^ + u^))- m'^{6Q^{i + uf+ 
Q^{5P + 3Pu + 3iu'^ + 5u^) - 2{P + hPu + ISPu^ + biu^ + u^)))). 
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2 4 6 8 10 12 14 

(GeV) 

Figure 1: The muon transverse momentum distribution da/dpi^ of the deep inelastic beauty 
production at HERA in the kinematic range > 2 GeV^ 0.05 <y < 0.7, p^^*^'"* > 6 GeV, 
-2 < ryj*^* < 2.5 and > 2 GeV, -0.9 < r/^ < 1.3 or > 2 GeV, -1.6 < < -0.9. The 
sohd, dashed, dash-dotted, dotted and short dash-dotted curves correspond to the predictions 
obtained with the J2003 set 1 — 3, KMR and KMS unintegrated gluon densities, respectively. 
The experimental data are from ZEUS [8]. 
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4 6 8 10 12 
pif (GeV) 

Figure 2: The muon transverse momentum distribution da/dpi^ of the deep inelastic beauty 
production at HERA in the kinematic range 2 < < 100 GeV^, Q.K y < 0.7, p^^^''''^ > 
6 GeV, Ir/j''*! < 2, > 2 GeV and -0.75 < 77^ < 1.15. Notations of all curves are the same 
as in Fig. 1. The experimental data are from HI [9]. 



17 








J L 



-1 







11 



Figure 3: The muon pseudo-rapidity distribution da/drj^ of the deep inelastic beauty pro- 
duction at HERA in the kinematic range Q"^ > 2 GeV^, 0.05 <y < 0.7, pJ^*^'*^^* > 6 GeV, 
-2 < Tyj^* < 2.5 and > 2 GeV, -0.9 < r]^' < 1.3 or p'^ > 2 GeV, -1.6 < 77'* < -0.9. 
Notations of all curves are the same as in Fig. 1. The experimental data are from ZEUS [8]. 
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Figure 4: The muon pseudo-rapidity distribution da/dr]^ of the deep inelastic beauty produc- 
tion at HERA in the kinematic range 2 < < iqo GeV^, 0.1 < y < 0.7, y^*^'^^* > 6 GeV, 
l^jctj < 2, > 2 GeV and —0.75 < rj^ < 1.15. Notations of all curves are the same as in 
Fig. 1. The experimental data are from HI [9]. 
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10 20 30 

^et Breit ^q^^j 

Figure 5: The jet transverse momentum distribution da / dj)^^^^^^^ of the deep inelastic beauty 
production at HERA in the kinematic range > 2 GeV^, 0.05 <y < 0.7, p^^*^"""'* > 6 GeV, 
-2 < Tyj^* < 2.5 and > 2 GeV, -0.9 < 77^ < 1.3 or > 2 GeV, -1.6 < 77'^ < -0.9. 
Notations of all curves are the same as in Fig. 1. The experimental data are from ZEUS [8]. 
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^et Breit ^q^^^ 

Figure 6: The jet transverse momentum distribution da / dp^^^^'^'^^^ of the deep inelastic beauty 
production at HERA in the kinematic range 2 < < 100 GeV^, 0.1 < y < 0.7, p^^^''''^ > 
6 GeV, Ir^j*"*! < 2, > 2 GeV and -0.75 < 77^ < 1.15. Notations of all curves are the same 
as in Fig. 1. The experimental data are from HI [9]. 
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Figure 7: The distribution of the deep inelastic beauty production at HERA in the 
kinematic range Q'^ > 2 GeV^, 0.05 < y < 0.7, p^^ > 6 GeV, -2 < rf^^ < 2.5 and 
> 2 GeV, -0.9 <r]'' < 1.3 or p'^ > 2 GeV, -1.6 < r]^' < -0.9. Notations of all curves are 
the same as in Fig. 1. The experimental data are from ZEUS [8]. 
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Figure 8: The distribution of the deep inelastic beauty production at HERA in the 
kinematic range 2 < < iqo GeV^, 0.1 < y < 0.7, p^^*^'*'^* > 6 GeV, |r/j^*| < 2, pf^ > 2 GeV 
and —0.75 < rj^ < 1.15. Notations of all curves are the same as in Fig. 1. The experimental 
data are from HI [9]. 
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lOQioX 

Figure 9: The log^Q x distribution of the deep inelastic beauty production at HERA in the 
kinematic range > 2 GeV\ 0.05 < y < 0.7, pJ^*^"""'* > 6 GeV, -2 < rf''^ < 2.5 and 
p^>2 GeV, -0.9 < Ty'* < 1.3 or p'^ > 2 GeV, -1.6 < r|^' < -0.9. Notations of all curves are 
the same as in Fig. 1. The experimental data are from ZEUS [8]. 
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lOQioX 

Figure 10: The log^g^ distribution of the deep inelastic beauty production at HERA in the 
kinematic range 2 < < iqq QeV^, 0.1 < y < 0.7, > 6 GeV, |r;j<^*| < 2, > 2 GeV 

and —0.75 < rj^ < 1.15. Notations of all curves are the same as in Fig. 1. The experimental 
data are from HI [9]. 
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Figure 11: The structure function FKa:,^^) as a function of x for different values of Q^. 
Notations of all curves are the same as in Fig. 1. The experimental data are from HI [6, 7]. 
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